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A general relativistic model of accretion disks with coronae 
surrounding Kerr black holes 

Bei You 1 , Xinwu Cao 1 , and Ye-Fei Yuan 2 
ABSTRACT 

We calculate the structure of a standard accretion disk with corona surround- 
ing a massive Kerr black hole in general relativistic frame, in which the corona 
is assumed to be heated by the reconnection of the strongly buoyant magnetic 
fields generated in the cold accretion disk. The emergent spectra of the accretion 
disk-corona systems are calculated by using the relativistic ray-tracing method. 
We propose a new method to calculate the emergent Comptonized spectra from 
the coronae. The spectra of the disk-corona systems with a modified a-magnetic 
stress show that both the hard X-ray spectral index and the hard X-ray bolo- 
metric correction factor £boi/-^x,2-iokcv increase with the dimensionless mass 
accretion rate, which are qualitatively consistent with the observations of active 
galactic nuclei (AGNs). The fraction of the power dissipated in the corona de- 
creases with increasing black hole spin parameter a, which leads to lower electron 
temperatures of the coronas for rapidly spinning black holes. The X-ray emis- 
sion from the coronas surrounding rapidly spinning black holes becomes weak and 
soft. The ratio of the X-ray luminosity to the optical/UV luminosity increases 
with the viewing angle, while the spectral shape in the X-ray band is insensitive 
with the viewing angle. We find that the spectral index in the infrared waveband 
depends on the mass accretion rate and the black hole spin a, which deviates 
from f v oc z/ 1 / 3 expected by the standard thin disk model. 

Subject headings: accretion, accretion disks, black hole physics, magnetic fields, 
galaxies: active 
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Introduction 



Active galactic nuclei (AGNs) are powered by the accretion of gases on to their cen- 
tral massive black holes, and the emission of the accretion disks is responsible for the ob- 
served multi-band spectral energy distributions (SEDs). It is suggested that the optically 
thick and geometrically thin accretion disks are present in luminous AGNs, and the ob- 



thin accretion disks (e.g.. IShakura &: Sunvaev 



Sun fc Malkan 1989; Laor fc Netzer 



1989 



Laor 



1973 



1990 



Shields 


1978 


Malkan & Sargent 


1982 


Czernv et al. 


2011; 


Hu & Zhang 


2012) 



One feature in the spectra of such standard thin accretion disks surr ounding massive black 
holes is, f v oc z/ 1//3 , in the infrared waveband ( Koratkar &; Blaeall999l ) . As the accretion disk 
extends to the region very close to the black hole, the general relativistic th i n accretion disk 
models were developed in the previous works (e.g.. iNovikov fc Thornd 19731 ; iPage fc Thorne 
1974 ; lLaor fc Netzer! Il989t iLi et al.l l2005t lAbramowicz fc Fragile! 1201 ll ) , which were applied 
successfully to fit the observed SEDs o f AGNs and the values of black hole spin parameter 
can be derived for some sources (e.g., lLaorl Il990l ; ICzerny et al.l 1201 ll ) . The values of the 



the theoretical accretion disk spectra 



spin parameter were estimated by the comparison o 
with the observed spectra in some X-ray binaries (iZhang et al.l Il997t ISteiner et al.l 12009 



Gou et al.ll2011t iKulkarni et al.l 1201 It iMcClintock et al.ll201ll ). However, the temperature of 
the standard thin accretion disks is too low to produce the observed power-law spectra of 
AGNs in the hard X-ray waveband, and the accretion disk - corona model was therefore sug- 



gested (e.g.. iGaleev et al.lll979t lHaardt fc Maraschil Il99ll ; ISvensson fc Zdziarskil Il994l ). In 



this scenario, the power-law hard X-ray spectra of AGNs are most likely due to the inverse 
Compton scattering of soft photons on a population of hot electrons in the corona above the 
disk. The soft photons are emitted from the cold disk, and most of them pass through the 
optically thin corona without being scattered, which emerge as the observed optical/UV con- 
tinuum emission of AGNs. A small fraction of the photons from the cold disk are Compton 



Svensson & Zdziarski 


1994; 


Kawaguchi et al. 


2001 



3 works (e.g., 


Haarc 


t & Maraschi 


Liu et al. 


20021 


Cao 


2009 


). Most 



gravitational energy is generated in the cold disk t hrough the turbulence, which is proba- 
bly produced by the magnetorotational instability fealbus fc Hawley 1991 ). A fraction of 



the energy generated in the cold disk should be transported to the corona. One of the most 
promising mechanisms is the energy of the strongly buoyant magnetic fields in the disk being 
tran sported vertically to heat the corona above the disk with the rec onnection of the fields 



e.g. iDi Matted 119981 iDi Matteo et al.l 119991 Merloni &; Fabian! l200lh . The temperature of 



the hot electrons in the corona is roughly around 10 9 K, w hich can successfully repro duce a 



power-law hard X-ray spectrum as observed in AGNs (e.g. iLiu et al.ll2003t ICaoll2009l ). 
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The emergent spectrum of an accretion disk surrounding a rotating black hole is in- 
fluenced by the relativistic effects due to the strong gravity fi e ld of the hole, which was 



extensively investigated in many works (e.g., IZhang et al.lll985t lLaor fc Netzerlll989l ; lLaor 



1991; 


Li et al. 


2005, 


2009) 



20091 ). The iron Ka fluorescence lines observed in some AGNs are sig- 
nificantly asymmetric, characterize d with a stee p blue wing, an extended red wing and the 
blueshifted emission line peak (see iMillerl 120071 . for a review, and the references therein). 
These features can be modeled as emission from the inner region of an accretion disk il- 
luminated by the external X-ray radiation with the Doppler boost effect and the general 
relativistic effects (frame -dragging, gravitational redshift and bending of the light) being 



properly considered (e.g., lLaorlll99lf ). Such effects should also be considered in the calcula- 



tions of the emergent continuu m spectrum of an accretion disk surrounding a rotat i ng black 



hole observ ed at infinity (e.g., lLaor fc Netzerl Il989l ; IZhang et al.l 119971 ; lLaorl Il99ll : iLi et al. 



20051 . 120091 ). The ray-tracing method is widely adopted to derive the photon trajectories 
bent by the strong gravity fiel d of the black hole in the calculations of the emergent spec 
tra emitted from the disk (e.g., 



1998; Fuerst fc Wu 2004; Yuan et al 



Chandrasekhar 



2009 



1983 



Rauch fc B 



Abramowicz &; Fragildl201ll ) . The values of black 



andford 1994; Cadez et al. 



hole spin parameter can be estimat ed by the compar ison of the observed spectra with the 



theoretical model calculations (e.g., I Yuan et al.ll2010l ). Besides the radiation from the disk, 



the emergent spectrum of the corona is also affected by these general relativistic effects, 
which is more complicated than the disk case due to the comple xity of the radiation transfe r 
of Comptonized photons in the geometrically thick corona (e.g.. ISchnittman fc Krolikl|2010l ). 



In this work, we calculate the structure of a thin accretion disk with corona surrounding 
a massive Kerr black hole in general relativistic frame, in which the corona is assumed to be 
heated by the reconnection of the magnetic fields generated by buoyancy instability in the 
cold accretion disk. The emergent spectra of the accretion disk-corona systems are calculated 
by using the relativistic ray-tracing method. The emergent Comptonized spectra from the 
coronae are calculated by dividing layers in the coronae. We summarize the calculations of 
the disk-corona structure and its emergent spectrum in Sections 2-4. The numerical results 
of the model calculations and the discussion are given in Sections 5 and 6. 



2. The accretion disk-corona structure 



In the Boyer-Lindquist coordinates (t,r,6,4>) with natural units G = c = 1, the Kerr 
metric is 

ds 2 = -e 2v dt 2 + e 2 ^(d<J) - udtf + e 2 ^dr 2 + e 2 ^d6 2 (1) 
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where 



,2k 



£A 



A 



2M bh r + a 2 



f sin 2 
E 



A' 6 



2^2 



2M hh ar 



r 2 + a 2 cos 2 9, £ 



,2\2 



a 2 A sin 2 9, 



Mbh is the mass of black hole, and a is the black hole spin parameter (— 1 < a < 1} 
event horizon of the black hole is 



M bh + (M, 



bh 



,2Nl/2 



(2) 
The 

(3) 



The Keplerian velocity of the matter moving around the black hole in circular orbits in the 
equatorial plane at radius R is given by 

n K = ± , hh 1/2 , (4) 

where the upper sign refers to direct orbits, i.e., corotating with the black hole spinning 
direction, and the lower sign refers to retrograde orbits, i.e., counterrotating with the black 
hole spinning direction. The innermost of a geometrically thin accretion disk is assumed to 
extend to the marginally stable circular orbit of the black hole, which is given by 



M bh {3 + Z 2 T [(3 - Z 1 )(3 + Z 1 + 2Z2)} 1 ' 2 } 



(5) 



where 



and 



Z l = l 



M 2 



1/3 



Mi 



3a 2 
Ml 



bh 



1/3 



1/2 



M 



bh 



1/3- 



In this work, we consider the relativistic accretion disk with corona, in which the 
cold disk is optically thick and geometrically thin. The circular motion of the matter in 
the accretion disk is assu med to be Keplerian, as that in the standard thin disk model 
(jShakura fc Sunyaevlll973l ). The formulation of the relativistic accretion disk model is sim- 
ilar to the standard thin disk mode l, in which some general relativistic correction factors 
A, B, C, D, and L, are included iNovikov fc Thorne! Il~973l : IPage fc Thornel Il974h . The 



values of all these quantiti es approach unity in the region far from the black hole (see 



Abramowicz &: Fragile] 1201 ll . for a review). 



The gravitational power of the matter dissipated in unit surface area of the accretion 
disk is 

3GM bh M L 
^ dissi 87T.R 3 BC 1 / 2 ' 



(6) 
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where M is the mass accretion rate of the disk. A fraction of the dissipated power is 
transported to the corona most probably by the magnetic fields generated in the disk. The 
fields are strongly b uoyant, and the c orona above the disk is heated with the reconnection of 
the magnetic fields fjDi Matteolll998l ). The soft photons emitted from the disk are Compton 
scattered by the hot electrons in the corona to high energy, which is the main cooling process 
in the corona. About half of the scattered photons are intercepted by the disk, part of whic h 
are reflected and the remainder are re-radiated as blackbody radiation (IZdziarski et al.lll999l ). 
Thus, the energy equation for the cold disk is 



1 



Qdissi-Qcc,+ o(W)Q 



4aT d 4 



(7) 



where Td is the interior temperature in the mid-plane of the disk, and r = r cs + is 
the optical depth in the vertical direction of the disk, and the reflection albedo / = 0.15 
i s adopted in all our calculations. The power transported from the disk to the corona is 
jDi Matted Il99~8h 

Qcoi 



B 2 



PmVp = —V p {$) 

where p m is the magnetic pressure in the disk, and v p is the vertically rising speed of the 
magnetic fields in the accretion disk. The rising speed v p is assumed to be proportional 
to internal Alfven velocity, i.e., v p = bv^, in which b is of the order of unity for extremely 
evacuated magnetic loops. The strength of the magnetic fields in the disk is crucial in our 
present investigation, which determines the ratio of power dissipated in the corona to that 
in the disk. This ratio can be constrained by the observed ratio of the X-ray to bolometric 
luminosities, which seems to support that the values of b should not deviate significantly 
from unity in most AGNs, though a low value of b (e.g. b < 0-1) cannot be ruled out 



in few individual AGNs with extremely low Lx/L^ Q \ (see IVasudevan fc Fabianl 120071 ; ICao 
20091 ). The detailed physics for generating magnetic fields in the accretion disks is still 
quite unclear, and we assume the magnetic st ress tensor to be related to the pressure of 
the d i sk, as done in many previous works ( e. g.. ISakimoto fc Coronitilll981t IStella fc Rosner 
19841 ; iTaam fc Lir]|l984l ; iffirose et al.l [2009h . ICaoT (]2009[ ) constructed accretion disk-corona 
models with different magnetic stress tensors, and found that the model calculations can 
qualitatively reproduce the observational features that both the hard X-ray spectral index 
and the hard X-ray bolometric correction factor £boi/-^x,2-iokev increase with the Eddington 



ratio, if r, 



rip 



a 



y/PgasPtot (Ptot = Pgas + Prad is adopted). We use this magnetic stress 



tensor in most of our calculations in this work. 



The angular momentum equation for the gas in the disk is (INovikov &; Thorndll973l ) 

(9) 



4irH d T rv = M 



GM hh F 



R3 D 
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where is the scale height of the accretion disk and the factor F is first introduced in 
Riffert fc Heroldl (1l995h in the form of a function D. In the vertical direction of the accretion 
disk, the hydrostatic equilibrium is assumed. The vertical pressure gradient being balanced 
with the vertical component of the gravity leads to (IRiffert fc Heroldlll995l ; iMcClintock et al. 
2006h 

dp pGM hh z J 
dz 



R 3 C 

which is valid for the geometrically thin disk, and J = 1 — 
of the disk is given by 

R 3 C 



where the sound speed c s 



4a I 3af 

r 3/2 T r 2 



(10) 

Thus the scale height 



Hd ^\GM hh J 

[(Ptot+Pm)/ P } 1/2 . 



1/2 



The continuity equation of the disk is 



M = -ATrRH d pV r D 1/2 , 



(12) 



where p(R) is the mean density of the disk, V r (R) is the radial velocity of the accretion flow 
at radius R, and the mass accreted in the corona is neglected. 



The state equation of the gas in the disk is 



Ptot — Pgas + Prad 



pkT d 
pm p 



(13) 



where ao is the radiation constant, the mean atomic mass p = + 1/Vc) l - We assume 
the plasma to consist of 3/4 hydrogen and 1/4 helium, i.e., \i\ = 1.23 and p e = 1.14. 

The ions and electrons in the corona are heated by the reconnection of the magnetic fields 
rising from the disk, of which the cooling is dominated by the inverse Compton scattering of 
the soft photons from the disk, the synchrotron radiation, and th e bremsstrahlung radiation. 
The energy equation of the two-temperature corona is given by (IDi Matted Il998[ ) 



Qcor Qc 



6Q- C 



(14) 



where F~ r = f S y n +^brem+^comp i s the cooling rate in unit surface area of the corona, and Q]? OT 
is the energy tran sfer rate from the ions to the electrons in the corona via Coulomb collisions, 
which is given by IStepney fc Guilbertl (119831 ) . We ado pt the fraction of the energy d i rectly 
heati ng the electrons S = 0.5 in all our calculations (jBisnovatyi-Kogan Sz Lovelace! 119971 . 



20001 ) . The cooling terms F and F hrcm are the functions of the magnetic field strength, 



the number density and temperature of the electrons in the corona, which are taken from 
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Narayan fc Yil (jl995[ ). The cooling rate -F Comp of the Compton sc a tterin g; in the corona is 
calculated by using the approach suggested by ICoppi fc Blandfordl (119901 ) . 



The state equation of the gas in the optically thin corona is 

Pcor kT\ Pcor kT e 



Pc 



/i e m p 



+ Pc 



(15) 



where T ; and T c are the temperatures of the ions and electrons in the two-temperature 
corona, and the magnetic pressure p cor ,m = B^ or /Sn. Given uncertainty of the magnetic 
fields in the corona, we sim ply assum e the magnetic pressure to be equipartition with the 
gas pressure in the corona (ICaoll2009i ). The cooling of the corona is always dominated by 



the inverse Compton scattering, and the strength of the fields in the corona mainly affects 
the synchrotron radiation. Our final results are almost independent of the value of the field 
strength adopted, except for the spectra in radio wavebands. The scale height of the corona 
H cor can be estimated based on the assumption of static hydrodynamical equilibrium in the 
vertical direction. 



As done in the work of ICaol (120091 ). we need to specify the ion temperature of the 
corona in the calculations of the disk-corona struct ure. It is found that the ion tempera- 
ture T{ has little influence on the X-ray spectra, and iLiu et al.l (120031 ) pointed out that the 
temperature of the ions in the corona is always in the range of ~ 0.2-0.3 virial tempera- 
ture (T V j r = G M\ ) \ i m v / kR) . In this work, the definition of the virial temperature is slightly 
different, 



3 -kT- 
2 vir 



-m p QlR 2 , 



(16) 



where fix is the Keplerian velocity of the gas in the disk at radius R. In this work, we 
ad opt T; = . 9T v j r in all our calculations, which is almost equivalent to the value suggested 
bv lLiu et all boO$ ). 



It is still quite unclear on the outer radius of the accretion disk. One conventional 
estimate of the outer radius of the disk is base d on the assumpt ion that the disk is truncated 
where it becomes gravitational unstable fe.g. . ICollin fc Zahnll2 



disk is estimated with the Toomre parameter ( Toomrel 

^kc s 



1964 



20081). The outer radius of the 



Goldreich fc Lynden-Belllll965l ) 



Q 



Toomre 



2nGpH d 



1. 



;i7) 



Given the values of the black hole mass Mbh, the spin parameter a, the viscosity pa- 
rameter a, and the mass accretion rate M, the disk structure can be derived by solving 
Equations ([6])- ([8]) numerically. The structure of the corona is then available based on the 
derived disk structure by solving Equations (fll| - (fT6l) . 
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3. Ray tracing method 



The trajectory of the photons emitted from the inner region of the accretion disk-corona 
system is bent by the strong gravity of the black hole. The observed emergent spectrum of 
such an accretion disk-corona system surrounding a Kerr black hole can be calculated with 
the ray tracing method. We summarize the method briefly in this section. 



In the Kerr spacet i me, t here are four constants for the motion of a photon (jCarterlll968 



Abramowicz &: Fragile! 1201 ll ). three of which can determine the trajectory of the photon, 
namely the energy at infinity, 

E = - Pt , (18) 



the effective angular momentum, 



and the Carter constant, 



(19) 



Q = p 2 e - a 2 E 2 cos 2 9 + L\ cot 2 9, (20) 

where p M is the four-momentum of the ph oton. Combi ning Equations O18p - (|20p and the 
expressions of a photon's four- momentum in lCarterl (119681 ). we derive the following equations 
describing the trajectory of a photon in the Kerr spacetime: 

, dr 



da 

E— 
da 

E— 
da 

E— 
da 



L 7 2arE 



L z a 2 



sin 2 9 A 
—a(aE sin 2 9 — L z 



a 2 )[E(r 2 + a 2 ) - L z a] 



where 



R(r) 
&(9) 



[ E {r 2 + a 2 ) - L z a] 2 - A[{L Z - aE) 2 + Q], 



Q _ (_ a 2£2 + 



T 2 
sin 2 ^' 



cos 2 9. 



(21) 
(22) 
(23) 
(24) 

(25) 
(26) 



Here a is the affine parameter along the trajectory. The signs in the above equations must 
be the same as the signs of dr, d9 respectively, and they change at the turning point in r 
or 9, i.e., R(r) = or Q(9) = 0. A s the trajectory of a photon is independent of its energy 
E, ICunningham fc Bardeenl (119731 ) introduce two dimensionless parameters being conserved 
along the trajectory, 

Q 



E' " ~ E 2 



(27) 
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and define another affine parameter dP = EjUda. The equations describing the trajectory 
of the photon can be re-written as 



dr 

de 

d<j> 
dt 



±dR(r)dP, 



±\/Q(9)dP, 

A 2ar — Aa 21 



sin 2 9 



dP, 



jr 2 + a 2 ) 2 
A 



- a 2 sin 2 9 - 



2arX 



dP, 



(28) 
(29) 
(30) 

(31) 



where R(r) = R(r)/E 2 , 0(6*) = Q(9)/E 2 . Integrating the null geodesic equations of the 
photon along the path from the initial location (r em , 9 em , em ) to the observer at infinity 
with position (r obs , obs , o bs), the trajectory can be determined by the following integral 
equations 

'"obs ^obs 

dr f d6 



± 



± 



R(r) 



P. 



9(0) 



and 



^obs 



± 



sin 2 fl 



+ 



2ar — a 2 X 



A 



dP, 



(32) 



(33) 



where P em , Pobs are functions of 9 or r, corresponding to the upper limit value of the 9- 
integral taken as 9 cm , 9 ohs respectively, and the signs in the above integrals must be the same 
as those of dr and d9 to guarantee the integral P is always positive and increasing along 
a photon's trajectory In the calculations, we set r D b s = oo, and o h s = 0. The detailed 



Cadez et al. 


1998; 


Li et al. 


2005) 



The image of the accretion disk-corona observed at an angle 9 Q ^ S with respect to the 
axis of the black hole at infinity can be described by the two impact parameters, 



a 



rp 



(</>) 



P 



(t) 



A 



sin 9, 



obs 



and 







rp 



(0) 



V 



(q + a 2 cos 2 # bs — A 2 cot 2 9, 



obs 



1 1/2 



Pobsj 



(34) 



(35) 



where a is the apparent displacement of the image perpendicular to the projected axis of 
the black hole, and (3 is the apparent parallel displacement. The vector is the four- 
momentum in the local nonrotating frame(LNRF). The observed image of the accretion 
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disk-corona can be calculated with the ray tracing method described above, provided that 
the local spectrum of the disk-corona is known. For a given point (a, (3) in the image seen 
by the observer at (oo,# o b s ,0), the two constants of motion A and q can be derived with 
Equations 034p and (I35p . Then the trajectory integral P, the radius r em , and the azimuthal 
angle em , can be calculated by using Equations (|32|) and (|33|) . 



4. Radiation transfer and emergent spectrum 

In this work, the photons cannot move through the equatorial plane because of the 
optically thick accretion disk, and therefore we neglect such trajectories. For the photons 
passing through the corona, the optical depth for the scattering of photons in the corona in 
the curved space can be calculated with the ray tracing method. Along the trajectory of the 
photon in the corona, the differential optical depth is 

dr = aTn e dl, (36) 

where or is the Thompson cross-section, n e is the electrons number density in the corona at 
(R, z), and 

dl = e v Y,dP = ± e Tj dO. (37) 
Integrating along the photon trajectory of the photon in the corona, we have 

r = ± T ^£3 de, (38) 
0(9) 



where d6 = ±^Q(9)dP. 

As the accretion disk is geometrically thin, we approximate 8 cm = n/2 in calculating 
the emergent spectrum of the disk. The radiative flux from the accretion disk observed at 
infinity is given by 

dad(3 

IT 2 

where B Uobs is the intensity at observed frequency z/ Q b s , dtt = dad(3/D 2 is the solid angle 
subtended by the image of the disk. 



FKbs) = / e-B Uohs dn = I e-B Uobs —f, (39) 



The gravitational redshift effect and Dopple r redshift can be included by using the 
Lorentz invariant I/u 3 ( IRybicki Sz Lightmanlll979l ). so the emergent spectrum can be calcu- 
lated with 



' 9 3 e- T B Uc ^ (40) 
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where B Ve is the the intensity of the photons from the disk at emitting frequency v e . Inte- 
grating over the disk image, the emergent spectrum of the disk is available with 



WKbs) = ^D 2 F(u ohs ) = An / g 6 e T B Uc dad/3 



(41) 



The redshift factor is defined by ([Cunningham fc Bardeenlll973l ) 



^obs PfjV' lobs 



(42) 



where is the four- momentum of the photon (jCarterlll968l ). 



Pft = {Pt,Pr,Pe,P<t>) = 



R(r) 



A 



±yjQ(9),X)E. 



(43) 



The four-velocity of the fluid in the locally non-rotating frame (LNRF) u^ hs and the local 
rest frame (LRF) u^ m are given by 



obs 



1,0,0,0), 



and 



(44) 



(45) 



respectively ( lYuan et al.ll2009l ) 



The corona is assumed to rotate with the disk at Keplerian velocity. The calculations 
of the emergent spectrum of the corona is slightly different from that of the accretion disk, 
because the redshift factor g and the trajectory of a photon depends on the location of the 
photon emitted (r em , em , z em ) in the corona. To calculate the redshift factor g and the 
trajectory of such a photon, we numerically solve Equations (|32|) and (133]) to derive the 
impact parameters a and /3. We divide the corona into small volumes with dhde where 
de is element area and dh is the thickness of a layer, and the emergent spectrum of each 
volume is available by calculating the radiation transfer with the ray tracing method, if the 
emissivity of the corona is available. For simplicity, we assume the radiation of the corona 
to be isotropic and the emissivity is homogeneous in the vertical direction of the corona. 



The synchrotron and bremsstrahlung emissivities are taken from lNarayan Sz Yil (119951 ). 
The spectrum of synchrotron and bremsstrahlung emission from the corona is available by 
summing up the contribution from these small volumes, 



7 syn,brem 



47T 



syn + £ brem )dad ^ 



(46) 
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where the effects of the absorption and the scattering in the corona are considered as de- 
scribed at the beginning of this section. 

In the calculati ons of the Comptonization in the corona, we assume the corona to be 



a parallel plane (see ICaol 12009k for the details). The mean probability of the soft photons 



injected from the cold disk experiencing first-order scattering in the corona is 

i 

Pi = 2 f(l- e~ To/cose ) cos9dcos9, (47) 



o 



where the constant specific intensity of the soft photons from the disk is assumed, 9 is the 
angle of the motion of the soft photons with respect to the vertical direction of the disk, 
and r = a T n c H cor is the optical depth of the corona for electron scattering in the vertical 
direction. We can calcul ate the first-order Compto nized spectra of the corona F^° mp using 



the method suggested by ICoppi fc Blandfordl (jl990l ) with the probability given by Equation 



(j4T|) . if the density and the temperature of electrons in the corona, and the incident spectrum 
of the cold disk are known. As we have assumed the scattered photons to be isotropic, 
the emissivity of the first-order Compton scattered photons in the layer with the height 
between h and h + dh can be approximated as e^° mp (i? d ) = F^ omp dh/ H COT . The first- 
order Compton scattered photons will pass though the corona before arriving the observer. 
Using the above described ray tracing approach, we can calculate the emergent first-order 
Comptonized spectrum of the corona. 

The mean probability for these first-order scattered photons experiencing the second- 
order scattering can be estimated as 

i i 

P 2 = I J df y [1 - e^ To/cose + 1 - e- (w)To/cos9 ]dcos#, (48) 



where £ = z/H cor , and the first-order scattered photons are assumed to be radiated isotrop- 
ically. The two terms in the integration are for downward- moving and upward-moving pho- 
tons emitted at height z, respectively. The probability of the scattered photons experiencing 
next higher order scattering approximates to P 2 , so we simply adopt P n = Pi for n > 2. 
Thus, the nth-order Comptonized spectra -F^° mp and emissivity e[;° mp can be derived with 
the same method for calculating the first-order Comptonization, in which the (n — l)th-order 
Comptonized spectra are used instead of the incident spectrum from the cold disk used for 
the calculation of the first-order Comptonization. In the same way, the emergent spectra 
of high order Comptonized spectra are derived, and the total Comptonized spectrum of the 
corona is available by summing up the contribution from the whole corona, 

iccnp = 47T J ^( £ C°mp + £ Comp + + £ Corn P)dad ^ (4g) 
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5. Results 

As described in Section 2, we can calculate the structure of an accretion disk with corona 
surrounding a rotating Kerr black hole, if the black hole mass Mbh, the spin parameter a, 
and the dimensionless accretion rate rh are specified. Based on the derived structure of the 
accretion disk-corona system, the emergent spectrum observed at infinity with an inclination 
angle o b s is calculated with the ray tracing method described in Sections 3 and 4, in which all 
the general relativistic effects are considered. The black hole mass M bh = 1O 8 M is adopted 
in most of our calculations for AGNs. The dimensionless mass accretion rate is defined as 
rh = M/M Edd (M Edd = L Edd /0.1c 2 ) independent of black hole spin. 

We plot the ratios of the power dissipated in the corona to the total released gravitational 
power L cor / Lboi as functions of the accretion rate rh with different values of the black hole spin 
parameter a in Figured] It is found that the dependence of the ratio L cor /Lboi with rh is quite 
different for the cases w ith d i fferent magnetic stress tensors, which are qualitatively consistent 



with those obtained in ICaol (120091 ) for the Newtonian accretion disk corona model. We find 
that the ratio L cor /Lboi decreases with increasing spin parameter a with either r rip = aP gas 
or r Tip = a^/P gas P tot if all other parameters are fixed, which means that the efficiency of 
the energy transportation from the disk to the corona is low for a rapidly spinning black 
hole independent of the detailed magnetic stress tenor adopted. However the model with 
r rip = aP tot appears to show the opposite trend that L cor /L bo i reaches large values for a 
rapidly spinning black hole. 

In Figure [21 we plot the interior temperatures of the accretion disks as functions of 
radius R with different values of the accretion rate rh and black hole spin parameter a, and 
the stress r Tip = a *J P gas -Ptot is adopted in the calculations. It is not surprising that the 
temperature of the disk increases with the mass accretion rate rh. For a rapidly spinning 
black hole, the inner edge of the disk extends to the region close to the black hole, and the 
temperature of the gas in the inner edge of the disk is significantly higher than that for a 
slowly spinning (or non-spinning) black hole. 

The structures of the corona calculated with r Iip = cty/P gas P tot , i.e., the temperatures 
of the ions and electrons, and the electron scattering optically depth of the corona in the 
vertical direction, are plotted in Figure [3j The electron temperature decreases with increasing 
mass accretion rate rh and the black hole spin parameter a. The optical depth for electron 
scattering is in the range of ~ 0.1 — 0.3 for different values of rh and a adopted. 

The emergent spectra of the accretion disk-corona systems with r rip = «-y/P gas P tot ob- 
served at infinity with different disk model parameters and viewing angles are plotted in 
Figure HI In each panel, one of these parameters is fixed at a certain value and the other 
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ones are taken as free parameters in the calculations of the emergent spectra. The corona 
is optically thin and geometrically thick, and therefore its radiation is almost isotropic. The 
radiation of the disk is dominantly in the optical/UV bands, the emergent spectrum of which 
depends on the viewing angle sensitively due to its slab-like geometry. The photon spectral 
indices, and the ratios of the bolometric luminosity to the X-ray luminosity in 2 — 10 keV 
as functions of the accretion rate rh for different values of a and 6* b s are plotted in Figure 
[5j The hard X-ray spectral index T increases with the accretion rate rh and black hole spin 
a, while it is insensitive with inclination # b s . The ratio of the bolometric luminosity to the 
X-ray luminosity increases with mass accretion rate, which is higher for a rapidly spinning 
black hole. 

The calculated emergent spectrum is usually insensitive to the outer radius of the ac- 
cretion disk, except the spectral shape in the infrared/optical wavebands. The outer radius 
of the accretion disk can be estimated by assuming the disk to be truncated where the self- 
gravity of the disk dominates over the central gravity of the black hole. The outer radii of 
the accretion disks with different accretion rates rh and black hole masses Mbh are plotted 
in Figure O The results show that the accretion disk surround a less massive black hole 
extends to a larger radius. The spectral shape of the accretion disk in the infrared/optical 
wavebands also varies with the outer radius of the disk -R ut- We plot the emergent spectra 
of the disk-corona systems with r Ttp = ay/P gas P t ot in the near- infrared/optical wavebands 
with different values of the accretion rate rh and black hole spin parameter a in Figure We 
find that the spectral shapes of the accretion disk-corona systems devia te significantly from 
f v oc z/ 1 / 3 , which is predicted by the standard a cc retion disk model (jshakura h Sunyaev 



19731 : iKoratkar & Blaesl Il999l : iHubenv et all 120001 ). iDavis et all (120071 ) measured UV spec 



tral slopes (s, where /„ oc i> s ) for several thousand quasars from the Sloan Digital Sky Survey 
(SDSS), and the distributions of the slope s span a wide range rather than s = 1/3 predicted 
by the standard accretion model. 



6. Discussion 



In this work, we calculate the structure of an accretion disk with corona around a Kerr 
black hole, and the emergent spectrum of the system is derived with the ray tracing approach, 
in which the general relativistic effects have been considered. The ions and electrons in the 
corona are assumed to be heated by magn etic power with the reconnection of the buoyant 
magnetic fields transported from the disk (IDi Matted Il998l ). The strength of the magnetic 
fields can be estimated from gas/radiation pressure of the disk with assumed magnetic stress 
tensor. We examine the relations of the ratio L cor /Lboi with the mass accretion rate rh in 
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the models with different magnetic stress tensors (see Figured]). In order to understand the 
results, we analyze the situation based on the standard thin disk mod el. The gravitational 
powe r dissipated in unit surface area of the accretion disk is given by (IShakura fc Sunyaev 
1973h 



Ha 



- r ripH d Qi, 



n T ripC s 



Ptot + Pn 
P 



1/2 



(50) 



where W T(p = f r Tip dh, p to t = Pgas + Pr&d- The power dissipated in the corona is 



PmVr, = bVAPv 



bp m B 
(4vrp) 1 /2 



bp* 



2gn 
P 



1/2 



(51) 



Combining Equations floTJ]) and (l5Tp . we find that the ratio L cor /L ho \ oc p m /(p m + Ptot)- If the 

aptot is adopted, the ratio L cov /L\ M \ remains constant independent 
'up = Pm = apgas, or T Iip = OLyjPgasPtoti the ratio 



stress tensor T T(p = p n 
of accretion rate rh. 



For the cases of r, 

-3/2 



Lcor/Lboi oc or T^ 0/ * respectively, as p gas ~ T d , p rad ~ and the radiation pressure is 
dominant in the inner regions of the accretion disk. As the disk temperature increases with 
mass accretion rate rh, the ratio L COT / L^ Q \ therefore decreases with increasing rh. The results 
ob tained in th is work for spinning black holes are qualitatively consistent with those derived 



in 



Caol fl2009h . 



As the value of black hole spin parameter a increases, the disk extends close to the black 
hole, and more gravitational power is released in the inner region of the accretion disk. The 
temperature of the inner disk region T d increases with the spin parameter a, which increases 
the energy density of the soft photons emitted from the disk. This leads to the corona cooling 
efficiently, because the cooling of the corona is dominated by the Compton scattering of the 
soft photons from the disk. The electron temperature of the corona in the inner region of 
the disk is therefore lower for rapidly spinning black holes (see Figure [3]). 

The resulted spectra of the disk-corona systems with magnetic stress tensor r rtp = 
a^PgasPtot show that both the hard X-ray spectral index and the hard X-ray bolomet- 
ric correction factor Lboi/£x,2-iokeV increase with the Eddington ra t io, which are qualita- 



tively consistent with the observations of AGNs ( Wang et al.l 20041 ; IShemmer et al.l 12006 



200 



The observed ratio of 



Vasudevan fc Fabian! 120071 ) , and X-ray binaries (jWu fc Gu 
L np t / Lhni may provide additional information of the disk-corona sys tems fjDavis fc Laor 
201 lh . The stress tensor r rtp = a^/PgasPtot wa s initially suggested by iTaam fc Linl (119841 ) 



based on the idea that the viscosity is proportional to the gas pressure, but the size of tur- 
bulence should be limited by the disk scale height determined by the total pressure. This is 
also supported by the analysis on the local dynamical instabilities in magnetized, radiation- 



- 16 - 



pressure-supported accretion disks (IBlaes fc Socratesl 1200 It iMerloni fc Fabianl 120021). How- 
ever, recent radiative MHD simulations seem to prefer the stress tensor r rip = ap to t (IHirose et al 
20091 ). Our model calculations can be carried out when the form of magnetic stress tenor is 
specified. For simplicity, in most cases of this work, we calculate the structure and emergent 
spectra of accretion disk-corona systems using the magnetic stress tensor r rip 



a 



y/PgasPtot- 



The radiation from the corona is almost isotropic due to its optically thin and geomet- 
rically thick structure, while the radiation from the disk is anisotropic. It is therefore found 
that the hard X-ray spectral index T is insensitive with the viewing angle # b s - The X-ray 
spectra are mainly dominated by the inverse Compton scattering, which mainly depends on 
the temperature of the electrons. Therefore the trend of the photon spectral indices T and 
the X-ray correction factors -^ bol /L Xi 2-iokeV can be explained by the decrease of the tempera- 
ture of the electrons with increasing mass accretion rate m, or the black hole spin parameter 
a (see Figure [51 and also Figure |3]). The photon index T of the hard X-ray spectrum becomes 
large when the black hole is rapidly spinning if all other parameters are fixed, which is due to 
the temperature of the electrons decreases with increasing value of a. The correction factor 
£boi/£x,2-iokev derived by the model calculations with r rip = a^/v^Vtnt are significantly 



large r than the observed values when m ~ 0.1 (e.g., IWang et al.ll2004t IVasudevan fc Fabian 
20071 ). One possibility is th at an advectiori dominated accretion flow (ADAF) is present in 
the inner region of the disk (INarayan fc Yilll994l. Il995f). The inner A DAF is X-ray luminous, 
which leads to a low value of Lboi/£x,2-iokeV ( jQuataert et al.lll999l ). This is also consistent 
with the fact that th e X-ray sp ectra of the AGNs are hard when accreting at low rates 



see 



Wang et al.l l2004j ; ICaol 120091 for the detailed discussion). An alternative resolution is 
the model with magnetic stress tenor r rtp = ap to t, however, it always predicts a constant 
£boi/£x,2-iokev independent of mass accretion rate fa. This needs further investigations in 
the future. 

In Figure [7J the spectral shapes of the disk-corona model with different values of accre- 
tion rates fa and black hole spin parameter a in the near- infrared/optical bands are plotted, 
which is dominantly contributed by the blackbody radiation from the disk. The spectral 
shapes are significantly different from f v oc z/ 1 / 3 , as expected in standard thin disk model. 
This is mainly caused by the independent temperature distributions in the disks surrounding 
massive black holes and the outer radii of the disks being different from those in standard 
thin disk model. In the standard thin disk model, the spectral shape f v oc z/ 1 / 3 in the 
near- infrared/optical wavebands is derived by assuming the temperature of the accretion 
disk T<j ~ i?~ 3 / 4 and the disk to extend to infinity. However, the temperature profiles of 
relativistic accretion disks deviate significantly from ~ R~ 3 ^ in the inner region of the 
disk (R < 100i?s) (see Figure [2]). The contribution to the spectra in IR/optical bands from 
the blackbody radiation of the outer region of the disks can also be quite important. The 
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outer radius of the disk is a function of black hole mass Mbh and accretion rate m (see 
Figure [6]). The observati onal slope of the infra red continua can be used to constrain the 
accretion disk model (see iKishimoto et al.ll2008l ). It is not suitable to compare the observed 
spectra only with, f u oc v 1 / 3 , as predicted by the standard accretion disk model. In this 
work, the local radiation from the disk is approximated as blackbody emission. A better 
approximatio n is to include a spectral hardening factor in calculating the local spectrum of 
the disk (e.g.. IShimura fc Takaharalll995l ; iHubeny et al.ll2000l ; iDavis et al.ll2007l ). with which 
the energy peak of the calculated disk spectrum will shift to a slightly higher frequency, while 
the spectral shape changes little. This would be important for detailed fitting the observed 
SED of AGNs, which is beyond the scope of this paper. 

It is assumed that zero torque at the inn er edge of the disk ( Shakura fc Sunvaevl 



which was doubted in some previous works flAgol fc Kroliki |2000| ; lHawley fc Krolik 



1973). 



200lh . 



The general relativistic magnetohydrodynamic (G RMHD) simulations provide an effectiv e 
tool to explore the torque in the plunging region ( iShafee et all 120081 ; iKulkarni et all l201ll ) . 
In the accretion disk with the thickness h/r ~ 0.05 — 0.1 the magnetic torque at the radius 
of the innermost stable circular orbit (ISCO) is only ~ 2% of the inward flux of angular 
momentum at this radius, w hich indicates that t he zero torque is re ally a good approximation 

( 120121 ) investigated the accretion 



for geometrically thin disks ( IShafee et al. I l2008h . iGong et al 



disk-corona model similar to that developed by ICaol (120091 ). in which a non-zero torque at 



the inner edge of the disk is assumed. As the detailed accretion physics in the plunging 
region is complicated and still quite unclear, we adopt the assumption of no torque in the 
inner edge of the accretion disk in all our calculations. The calculations in this work were 
carried out for massive black holes, which will be compared with the observed spectra of 
AGNs in our future work. The model calculations can also be applicable to X-ray binaries. 
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Fig. 1. — The ratios L cor /Lb i as functions of accretion rate rh with different magnetic stress 
tensors. The color lines represent different values of black hole spin parameter, i.e., a = 
(blue), 0.5 (green), and 0.998 (red). The different type lines correspond to different magnetic 
stress tensors, i.e., r rip = ap to t (dotted lines), r Tip = a^PgasPtot (solid lines), and r rip = ap g3 , s 
(dashed lines). 
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Fig. 2. — The temperatures in the mid-plane of the accretion disks as functions of radius 
R in units of R g = GM^/c 2 with different values of accretion rate rh and black hole spin 
parameter a. The colored lines represent the model results with different values of black hole 
spin parameter, i.e., a = (blue), 0.5 (green), and 0.998 (red). The solid lines represent the 
results calculated with rh = 0.01, while the dotted and dash lines are for rh = 0.1 and 0.5, 
respectively. The stress r T!fi = a^P gas P to t is adopted. 
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Fig. 3. — The upper panel: the temperatures of the ions in the coronae as functions of radius 
R with r rip = a y/Pg a s-Ptot- The middle panel: the temperatures of the hot electrons in the 
coronae. The lower panel: the optical depth for the Compton scattering of the electrons in 
the vertical direction of the coronae. The colored lines represent different values of black 
hole spin parameter, i.e., a = 0.998(red), 0.5 (green), and (blue). The solid lines represent 
the results calculated with rh = 0.01, while the dotted and dash lines are for m = 0.1 and 
0.5, respectively. 
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KHz) 

Fig. 4. — The spectra of the disk-corona systems with different values of rh, a, and # bs- The 
upper panel: the emergent spectra of the disk-corona systems surrounding spinning black 
holes observed at an angle 8 Q ^ S = 30° with respect to the disk axis (red: a = 0.998, green: 
a = 0.5, and blue: a = 0). The different type lines represent the spectra calculated with 
different values of accretion rate, i.e., rh = 0.5 (dash-dotted), 0.05 (dashed), and 0.01 (solid). 
The middle panel: the spectra of the disk-corona systems surrounding black holes spinning 
at a = 0.5. The different colored lines represent the emergent spectra observed at different 
inclination angles, i.e., 9 ohs = 60° (red), 30° (green), and 10° (blue). The different type lines 
represent the different values of accretion rate, i.e., rh = 0.5 (dash-dotted), 0.05 (dashed), 
and 0.01 (solid). The bottom panel: the spectra of the disk-corona systems calculated with a 
given accretion rate, rh = 0.1. The different colored lines represent the spectra with different 
values of black hole spin parameter, i.e., a = 0.998 (red), 0.5 (green), and (blue). The 
spectra observed at different angles are plotted with different type lines, i.e., # b s = 60° 
(dash-dotted), 30° (dashed), and 10° (solid). The stress r Tip = a^P gas P tot is adopted in the 
calculations. 
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Fig. 5. — The upper panel: the photon spectral indices as functions of accretion rate m for 
different values of black spin parameter a and viewing angle ^obs with r Tlp = a^/P gas P tot . The 
lower panel: the X-ray bolometric correction factors £t, i/£x,2-iokev as functions of accretion 
rate m for different model parameters. The colored lines represent the results calculated with 
a = 0.9 (red), 0.5 (green), and (blue), respectively. The different type lines represent the 
spectra observed at different angles, i.e., # b s = 60° (dash-dotted), 30° (dashed), 10° (dotted). 
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Fig. 6. — The outer radii of the accretion disks as functions of accretion rate m for different 
black hole mass M bh with r I(p = a *J P gas Ptot (blue: M hh = 1O 8 M , green: M hh = 5 x 1O 8 M , 
red: M bh = 1O 9 M ). 
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Fig. 7. — The spectral shapes of the disk-corona systems in the near IR/optical wavebands 
with T Tip = a<J '-Pgas-Ptot- The colored lines represent the model calculations with different 
values of black hole spin parameter, i.e., a = (blue), 0.5 (green), and 0.998 (red). The 
different type lines represent the results calculated with different values of accretion rate, i.e., 
m = O.Of (solid), O.f (dotted), and 0.5 (dash-dotted). The dashed line represents f v oc z/ 1 / 3 , 
expected by the standard thin disk model. 



